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Abstract 

The consequences of a recently introduced irreducible pionic effect in low 
energy nucleon-deuteron scattering are analyzed. Differential cross-sections, 
nucleon (vector) and deuteron (vector and tensor) analyzing powers, and 
four different polarization transfer coefficients have been considered. This 
3NF-like effect is generated by the pion-exchange diagram in presence of a 
two-nucleon correlation and is partially cancelled by meson-retardation con- 
tributions. Indications are provided that such type of effects are capable to 
selectively increase the vector (nucleon and deuteron) analyzing powers, while 
in the considered energy range they are almost negligible on the differential 
cross sections. These indications, observed with different realistic nucleon- 
nucleon interactions, provide additional evidences that such 3A^F-like effects 
have indeed the potential to solve the puzzle of the vector analyzing powers. 
Smaller but non negligible effects are observed for the other spin observables. 
In some cases, we find that the modifications introduced by such pionic ef- 
fects on these spin observables (other than the vector analyzing powers) are 
significant and interesting and could be observed by experiments. 
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I. INTRODUCTION 



Three- nucleoli forces as truly genuine interactions are not made by nature Rather, 
they are artifacts of theoreticians who decided to treat nuclear dynamics within a restricted 
Hilbert space. 

A sufficiently detailed picture for most low- and intermediate-energy phenomena in nu- 
clear dynamics should involve as constituents the nucleons, isobars and mesons (pions and 
heavier bosons). In doing so, all subnucleonic degrees of freedom, like quarks and gluons, 
are transformed into the notion of collective motion of hadronic clusters, which correspond 
to baryons and mesons. Such a system does not need the concept of two-, three-, and 
many-nucleon forces, since the mechanisms of meson emission and absorption, the ones that 
mediate the nuclear interactions, appear in the theory explicitly. 

For systems with more than two nucleons, not all the relevant hadronic degrees of freedom 
can be traced back explicitly in practical calculations, and one is faced with the problem 
to restrict the Hilbert space. This restriction generates three- (or many-) nucleon forces. 
And, most importantly, three-nucleon forces of different structure and of different physical 
interpretation arise depending on the different theoretical strategies which have been adopted 
to restrict the original hadronic system into the simplified, more tractable, Hilbert space of 
few-nucleon systems. 

A paradigmatic example that clarifies this situation is the one discussed in Ref. 
where two different strategies for describing the three-nucleon force are considered. Both 
approaches make the assumption that meson-exchange phenomena occur instantaneously 
and hence they both ignore any additional mesonic aspects brought in by the dynamics of 
the coupled meson-baryon system. The difference between the two considered approaches 
is that one of the two treats explicitly the non-nucleonic degree of freedom of the A isobar, 
thereby allowing the isobar to propagate over a finite time interval, while the other deals with 
nucleonic degrees of freedom only (the standard nuclear theory approach) which implies that 
one has to consider effective three-nucleon interactions explicitly, including those effective 
contributions for not allowing the isobar to propagate freely into the system. Following 
the first strategy, one ends up dealing with effective three-nucleon force {3NF) diagrams 
which are energy dependent due to the isobar propagation, and in addition one finds that, 
for consistency, also the two-nucleon interaction is modified in an energy dependent way 
because of the the presence of the third nucleon. The two-nucleon interaction in particular 
is modified in those contributions where an intermediate isobar propagation occurs, since the 
A is allowed to adjust to the nuclear medium by the dynamical approach. This secondary 
effect on the 2N force (sometimes called "dispersive 2N effect" ) counteracts the contribution 
due to the A-mediated 27r-exchange 3NF, and the physical interpretation of its origin is 
well understood in the approach with explicit A isobars. 

In contrast, the second strategy uses two- and three-nucleon interactions which are not 
allowed to adjust to the surrounding nuclear medium, because of the freezing of all the non 
nucleonic degrees of freedom. Conceptually, then it is very difficult to explain the origin of 
the dispersive 2A^ counter-terms, and the approach is amended to revive at least part of such 
contributions in an averaged way by introducing effective 3NF counterterms of short-range 
type. 

This paradigmatic example suggests that in the standard "few-nucleon" approach, with 
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its restricted Hilbert space consisting of nucleons only, some additional 3A^F-like diagram 
could be missing and the most viable way to reveal these new terms is to work-out the nuclear 
dynamics starting from an enlarged space where the system has to be treated dynamically 
also in presence of mesons, especially of the pionic degrees of freedom. If the origin of 
the 3NF effects is due to the restriction to the nucleonic Hilbert space, then we have to 
consider carefully how we theoretically reach this description, in order to pin-down the 3NF 
contributions. This is one lesson we must learn from the dynamical treatment of the A 
isobar in the three-nucleon system. 

The first non-nucleonic ingredient that plays a crucial role for nuclear systems is the pion, 
and it is natural to aim at a theoretical formulation that couples dynamically this meson 
to the nucleonic degrees of freedom. Such a formulation would then provide a combined 
description of nuclear systems at low- and intermediate energies. Pions are produced indeed 
and observed by experiments at intermediate energies and their dynamical role in these 
nuclear systems can hardly be ignored. 

With two-nucleon systems, various approaches to include the effects of a dynamical pion 
in a non-perturbative fashion have been developed and analyzed 0, and many theoretical 
problems like relativistic aspects, or the proper treatment of the nucleon renormalization 
effects, or finally the consistent treatment of meson-exchange diagrams related by different 
time orderings, have been extensively discussed ||^. Not all these problems have been solved 
in a systematic manner since these ttNN approaches aim to include just one dynamical pion 
in the theory, thus ruling out more complex situations with multipionic intermediate states. 
Stated in other words, these techniques allow to include one single dynamical pion in the 2N 
system, but then one has to face the conceptual problems arising because the states with 
more than one pion at the same time are ruled out from the dynamical equations. Still, in 
spite of these problems, the generalization of such type of approaches to the three-nucleon 
system allows to reveal possible new dynamical aspects of the pion that cannot be observed 
by freezing out from the very beginning all mesonic degrees of freedom, into instantaneous 
meson-exchange potentials. 

The dynamical equations for the coupled NN-tcNN problem have been generalized to 
the 3N system in Ref. 0. The method is based on the rigorous Yakubovski-Grassberger- 
Sandhas four-body theory, extended to the ttNNN system with an absorbable pion. 
Previously, this problem has been analyzed with few-body techniques for quite a few 
years P-[TT[|. The solution of the resulting set of 21x21 coupled equations represents a 
formidable task, and for practical reasons an approximation scheme to reduce the complex- 
ity of the original equations to an approximated but more tractable form has been devel- 
oped |[T^. The approximation scheme builds up the complete dynamical solution starting 
from the zeroth-order solution, given by the standard quantum-mechanical Alt-Grassberger- 
Sandhas (AGS) equation |T^. The strategy that emerges from this approximation scheme 
consists in treating the normal 2A^ correlations "exactly" via Faddeev-like methods, while 
the additional mesonic aspects which cannot be adequately described by a conventional 2A'^ 
potential have to be incorporated directly into the dynamical equation as corrections {i.e. 
through an underlying perturbative-iterative expansion). This is consistent with the find- 
ings of the approaches based on chiral perturbation theory (ChPT) which predict that 3NF 
effects are small 0. 

To the lowest order, the approximation scheme ends up with three different types of irre- 
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ducible 3NF diagrams which have to be incorporated in the dynamical equations according 
to the prescriptions given in Ref. [|l^]. Clearly, all three types of diagrams are related to 
specific aspects of the pion dynamics which cannot be incorporated in the description of a 
purely nucleonic system interacting through a pair-wise potential. 

Amongst the three types of diagrams, one can easily recognize the extensively studied 
2tt-3NF diagram. This diagram has led to various models of 3A^ potential amongst which we 
recall two historical representations, Tucson-Melbourne and Brasil |]2D|. These potentials 
use the nNN vertices and the off-shell extrapolated vrA^ amplitudes as inputs. It is important 
to observe that in the construction of these potentials certain diagrams have to be subtracted 
explicitly, to avoid double countings. The other two types of diagrams that emerge at 
the lowest order involve the intermediate formation of a 2N correlation while one pion 
is being exchanged, and finally the intermediate formation of a 3A^ correlation, during the 
exchange process. The last type of diagram has been discussed in the previous literature only 
occasionally ||2T|, and to our knowledge, its effect has not yet been estimated quantitatively in 
realistic situations, although its relevance might not be so important because the probability 
that a full 3A^ correlation is formed during the meson-exchange process is expected to be 
small. 

The diagram involving the 2N correlation during a pion-exchange process has been dis- 
cussed in Ref. |Ti], and the consequent appearance of an irreducible 3A^ operator with 



tensor-like structure has also been observed therein (see also Refs. |]T5|,|l6[). In the con- 
struction of this operator, certain classes of subdiagrams have to be subtracted, because 
of the presence of a cancellation effect which has been observed in the literature quite a 



few years ago |OJTH]. The connection between the nature of these cancellations and their 



implications to the ChPT approach has been discussed in Ref. [22|. A similar cancellation 



effect has been also observed - in leading order - in effective nuclear forces based on chiral 
lagrangians and constructed with the method of unitary transformation [^]. As has been 
discussed in Refs. [^^-[T^, however, such subtraction leads to a cancellation effect which 
is only partial if the original 3NF diagram has been derived within a dynamical approach 
which leads to an energy dependent 3A^F-like operator, and where the meson propagates 
for a sufficiently extended time to allow the intermediate formation of a correlated 2N pair. 
That might correspond to the inclusion of a series of non- vanishing higher-order terms in the 
chiral expansion. In the standard approach which uses instantaneous two- and three-nucleon 
interactions, the freezing out of the mesonic degrees of freedom does not contemplate the 
occurrence of such a 3NF effect. This situation presents interesting similarities with respect 
to that already encountered with the dispersive 2N counter-terms, which are properly taken 
into account only when the A degrees of freedom are treated explicitly in the dynamical 
equation. 

In Ref. [|l5l, the consequences of considering this type of one-pion-exchange 3NF diagram 
(0PE-3NF) in the 3A^ equation have been studied. It was shown that this diagram has the 
potential to modify considerably the vector analyzing powers Ay without affecting apprecia- 
bly the differential cross section, and could therefore solve the long standing puzzle of the 
vector analyzing powers in nucleon-deuteron scattering below 30 MeV. On the other hand, 
other, alternative explanations have been suggested in the literature. Effective three-nucleon 
potentials constructed with a combination of short-range and pion-range terms have been 
considered first under the point of view of the meson-exchange picture [E^, and later recon- 
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sidered under the framework of chiral perturbation theory p5[, since a non-neghgible role for 
these terms is predicted. A phenomenological spin-orbit three-body operator has also been 
introduced |^ in order to improve the vector analyzing powers. Another explanation p7 



advocates the Brown-Rho scaling-with-density hypothesis [^] in the three-nucleon system. 
Because of this scaling, the 2N potential has been modified specifically in the triplet P waves 
by reducing the scalar- and vector- meson masses to approximately 95% of their free-space 
value, thereby enhancing the spin-orbit term of the density-dependent 2N potential with 
respect to free space. This selective dependence of the Ay puzzle on the triplet P-waves 
of the 2N force suggested also that modern phase-shift analysis might have not yet been 
settled to the correct parameters for these states in the low-energy domain [^]. But it has 
also been argued that the required changes of the free-space 2N potential have to be 
exceedingly drastic for the one-pion-exchange contribution, which on the contrary is well 
established. 

In the present paper, we have extended the study of Ref. [T^ about the irreducible pionic 
effect implied by the 0PE-3NF diagram by considering such an effect with three different 2N 
potentials. In all cases we were able to demonstrate that this irreducible pionic effect has the 
potential to solve the puzzle for the low-energy vector analyzing powers with negligible effects 
for the differential cross sections and minor effects for the other spin observables. We have 
also widened the comparison with experimental data by including differential cross- sect ions, 
vector and tensor analyzing powers, and spin-transfer coefficients at various energies below 
20 MeV. 

The basic structure of the 0PE-3NF diagram is recalled in Sect.|I|. Comparison between 
theory and experiments is made in Sect. |IT1[ Conclusions are derived in Sect. 0. 



II. THEORY 



Following Refs. ||6| JT^Jl^JT5[| , we have incorporated directly into the Faddeev-AGS three- 
nucleon equation the irreducible effects generated by the one-pion-exchange mechanism in 
presence of a nucleon-nucleon correlation. The detailed expression has been discussed pre- 
viously in Ref. |M. It is reported here for convenience: 



^r(p,q,p',q;i^) 



(2vr) 



(1) 



(o-i ■ Q)(o-3 • Q)(ti ■ T3) + (0-2 ■ Q)(cr3 ■ Q)(t2 ■ T3) 



X 



ti2(p,p^;^-|; 

2m^ 



flMN^Q) 1 ti2(p,p';^- 17 -m, 

(27r)3 2m^ 



mt 



[<Ti ■ Q)(<T3 ■ Q)(ri ■ T3) + (0-2 ■ Q)(o-3 ■ Q)(r2 ■ T3) 



The momenta p, q represent respectively the Jacobi coordinates of the pair "12", and spec- 
tator "3" , while E is the 3iV energy. The pion-nucleon coupling constant is selected by the 
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underlying 2N potential that is considered; for instance, for the Paris and Bonn potentials, 
we have used the "traditional" value /^^^/(47r) = 0.078, while for the newer CD Bonn po- 



tential ( [^) we have consistently used the more recent determinations by the Nijmegen ||32 
and VPI group. Same considerations have been applied for the pion-nucleon form factor, 
since we have employed the same standard functions (and also with the same parameters) 
that have been employed at the level of the 2N potentials: 

fnNNiQ) = fnNN ^2 Q2 ' 



The transferred momentum Q = q' — q enters also in a;^ = ym^ + Q^. Uj denotes the 
subtracted t matrix between nucleons 1 and 2, defined according to the prescription 

ti2(p, p'; Z) = c{E) ti2(p, p'; Z) - vu{p, p') • (3) 
We have considered in addition another possible type of subtraction 

ti2(p, p'; Z) = c{E) ti2(p, p'; Z) - tu{p, p'; -A) , (4) 
where the subtraction parameter A has been fixed around 1.5 GeV. Other details can be 



found in Refs. |jT^,|T5|. The factor c{Z) is an adjustable parameter and serves to control 
the cancellation between the two terms. Ideally, this factor should be approximately one if 
the 2N t matrix could be reliably extrapolated off-shell down to Z ^ —160 MeV. However, 
the existing 2N potentials cannot guarantee the extrapolation at such negative energies. 
Moreover, additional approximations and simplifications entered in the determination of the 
expression for V^^ , as discussed in Ref. |]14|,|15[. For these reasons, we introduced c(Z) as 
adjustable parameter, with the constrain that at higher energies in nd scattering, this factor 
should move towards one, because the 2N t matrix entering in V^'^^ is then calculated at 
energies higher than ~ —160 MeV, that is, less far off shell. The use of the second type 
of subtraction, Eq. |^, has been introduced in the case of the Paris potential because this 
potential is not OBE-like (one-boson-exchange). In this case, the subtraction of the meson- 
exchange diagrams involved in the cancellation is therefore not feasible with Eq. ^. On 
the other hand Eq. ^ tends to enhance the cancellation effects with respect to Eq. ^, and 
therefore it is to be expected that the parameter governing the cancellation in Eq. ^ has to 
compensate this effect. 

The irreducible pionic effect described by Eq. ^ can be incorporated in the scattering 
equation in a convenient way if the 2N input potential is of finite rank. For a rank one (sep- 
arable) case, the 2N t matrix takes the expression t = \gi)T{gi\, and the (anti)symmetrized 
AGS equation can be reinterpreted as an effective two-body multichannel integral equation 
in one inter cluster momentum variable, 

Xii = Zn + ZiiT Xii, (5) 

with the driving term calculated as follows 

Zn = {9i\GoP\9i) + (gilGoVi'^'Golgi) . (6) 

The first contribution represents the standard AGS driving term, with Go and P being the 
free Green's function and the cyclic/anti-cyclic permutator, respectively, while the second 
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expression takes into account the effects of the irreducible 3A^ diagram discussed above. 
The procedure is consistent with the formahsm developed in Ref. 112] to include irreducible 
pionic effects as corrections in the Faddeev-AGS equation. Since the formalism is based 
on the systematic 4-body approach of Ref. it allows to include additional (high-order) 
classes of irreducible diagrams in subsequent steps. 

Neutron-deuteron scattering observables below 20 MeV have been calculated using the 
finite-rank representation of realistic nucleon-nucleon potentials, known as PEST, BBEST, 
and CDBEST potentials p4|-p6[]. These provide accurate representation of the nucleon- 



nucleon transition matrix for the Paris [^], Bonn B and CD-Bonn [^] potentials and 
are based upon the Ernst-Shakin-Thaler (EST) method for generating the finite-rank 
expressions of the transition matrices and/or potentials. Benchmark calculations for 
scattering and bound-state regime have demonstrated that with this method it is possible 
to solve accurately the Faddeev-AGS scattering equations, and to obtain results comparable 
(with errors of 1% or less) to those obtained from a direct solution of the 2-dimensional 
Faddeev equations (where the original 2A^ potential is used as input). The results shown 
in the next section have been calculated with the 2A^ potentials acting in the j < 2 states, 
and listed in Tab. |. For each 2A^ potential and state, the table also reports the rank of the 
separable expansion used. 

Finally, the calculations reported herein have been performed including all 3A^ states 
with total angular momenta up to J = 19/2, for both odd and even parities. 



III. RESULTS 

We first compare the results of our calculations with experimental data taken with inci- 
dent neutrons at 3 MeV. At this energy the process is below the deuteron break-up threshold. 
The experimental data shown in Fig. || are taken from Refs. |^ and |^3[ for the differen- 



tial cross-section and the neutron analyzing power, respectively. The figure compares the 
calculations obtained with the CD-Bonn, Bonn B, and Paris potential (dashed, dot-dashed, 
and dotted lines, respectively) with those obtained when adding consistently the 0PE-3NF 
effect as discussed in the previous section. By adjusting the parameter in Eq. ^ (or in Eq. ^ 
for the Paris potential), we obtain for each considered potential the required modifications 
for the proper reproduction of the analyzing power without ruining the description of the 
differential cross section. For the Bonn B and CD-Bonn potentials the actual values of the 
parameter have been reported in Tab. (For the Paris case the parameter has been set to 
0.385). 

For all three potentials the result of the calculations including the 0PE-3NF effect is 
contained in the thick solid line. For the differential cross-section, there is a very tiny 
effect or tendency to reduce the differential cross-section in both the forward and backward 
direction, however, such an effect can hardly be perceived in the figure. 

Then, we compare our results with other polarization observables, taken at comparable 
energies. To do so, we consider the complete set of deuteron analyzing power measurements 
in deuteron-proton scattering at 8 MeV, Ref. The measurements at this energy (for 



incident deuterons) compares to an equivalent energy of 4 MeV for the case of incident 
protons. However, we do not include Coulomb corrections in our calculations. As minimal 
Coulomb correction, we consider exclusively the additional loss of kinetic energy of the 
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proton while approaching the electric field of the deuteron . This loss amounts to about 
~ 0.7 MeV and hence we compare the experimental data of Ref. |44] with theoretical results 
obtained with incident neutrons of 3.3 MeV. We observe that the Coulomb slow-down effect 
has been observed experimentally for Ay and not for the deuteron analyzing powers, since 
there are no experimental neutron data here. However, a recent calculation [^] (with and 
without Coulomb force) suggests that the same effect holds also for these observables. 

The lower panel of Fig. ^ considers the deuteron vector analyzing power iTn. Again, 
we observe that the standard two-nucleon calculations without 3NF effects (thin lines) un- 
derpredict this observable, while the inclusion of the 0PE-3NF effect provides the necessary 
modifications (thick lines) suggested by the data. Conventions for the lines are the same as 
in the previous figure, namely, dashed, dot-dashed, and dotted lines, describe calculations 
with CD-Bonn, Bonn B, and Paris Potentials, respectively. We observe that the inclusion 
of the 0PE-3NF effect for the Bonn B potential provides the largest effect for iTn, while 
for the two other potentials the results with the 3NF effect are rather similar. This is at 
variance with respect to the case without 3NF effects, where the CD-Bonn potential alone 
provides a substantially higher iTu with respect to the Paris and Bonn B cases. Aside for 
the substantial increase of this observable provided by the 0PE-3NF effect, it is difficult to 
draw any additional conclusions by comparison with experimental data because of the pres- 
ence of the perturbation introduced by the Coulomb field, which at these energies modifies 
sensibly the shape of this observable not only in the most forward direction. 

The upper panel of Fig. ^ considers the deuteron tensor analyzing power, T22. Here the 
modification introduced by the 0PE-3NF are very small and cannot be perceived in the 
figure, with the exception of the Bonn B case, where a slight reduction of the dip peaked 
around 100° can be observed. A similar situation is suggested for the other two deuteron 
tensor analyzing powers, T20 and T21, shown respectively in the upper and lower panel of 
Fig. ^ The introduction of the 0PE-3NF effect introduces very slight modifications also 
in these two observables. Again we observe that with the inclusion of this 3NF effect the 
Paris and CD-Bonn results become more similar than they were before. Also with the Bonn 
B potential the modifications are small, although one can clearly perceive for T20 that the 
dip at 100° and the peak in the backward direction are slightly more pronounced. Same 
situation occurs for the dip at 80° for T21. Comparison with data shows that at these en- 
ergies and for these observables no definite conclusions can be drawn without an accurate 
inclusion of the Coulomb field, or without a comparison with accurate measurements in- 
volving neutron-deuteron scattering. The figure also suggests that with proper inclusion of 
Coulomb modifications in the theory and/or accurate neutron-deuteron measurements one 
could actually observe the phenomenological effects due to this 3A^F-like contribution. 

We have repeated the same analysis at 8.5 MeV. Fig. ^ exhibits that basically the same 
description found at 3 MeV holds also at this energy. The 3NF mechanism under scrutiny 
is capable to raise the neutron-deuteron Ay in order to match the experimental data without 
any sensible modification of the corresponding differential cross section. As a minor effect, 
also at this energy we can observe a reduction of the cross-section in the forward and 
backward directions, for both Bonn B and CD-Bonn potentials, but in both cases the effect 
is hardly perceptible in the figure. The cross-section data have been extracted from Ref. |^ 
(open circles) and Ref. [^] (open triangles) while the polarization data have been taken 
from Ref. 
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In Figs. 1^ and ^ we have compared "nc?" theoretical results with "j9(i" data for the 
deuteron analyzing powers. For the reasons explained above we compared results obtained 
at 8.3 MeV (for incident neutron energy) with data taken at 18 MeV (incident deuterons) 
5^ . The lower panel of Fig. ^ shows that when including the 3NF effect the increase of iTu 



is about of the right size for both potentials; however one observes also an inversion of the 
peaks, since, when including the 3NF effect, the Bonn B peak is higher than the CD-Bonn 
one, while without 3NF the CD-Bonn peak is higher. The same thing was occurring also 
at lower energy. 

For the deuteron tensor analyzing powers, T22 (upper panel of Fig. |^), T20 and T21 (upper 
and lower panels of Fig. ^, respectively), we observe that the overall shape of the angular 
distributions is better reproduced than at 4 MeV, indicating that the effects of Coulomb 
distortions are less important here (with the exception of the data in the forward direction), 
and that the modifications introduced by the 3NF effects are in general of minor entity with 
respect to those observed for the vector analyzing powers. A closer inspection, however, 
reveals that while the 0PE-3NF effects are very small for T22 for both CD-Bonn and Bonn 
B potentials, the situation is different for the other two deuteron tensor observables, where 
the effects of this 3NF diagram can indeed be observed and are about of the same size of 
the difference between the two potentials themselves. For T20 in particular, the 3NF effect 
provides a remarkable improvement of the dip at 110° in the case of the Bonn B potential, 
while for the case of the CD-Bonn potential the situation is basically unchanged. For T21 
again the 3NF effect remarkably improves the description in the Bonn B case (especially at 
the dip around 90°) , while for the case of the CD-Bonn potential the situation is reversed. 

Then, we have considered how the situation evolves at 12 MeV. Again we find that the 
introduction of these irreducible pionic effects are able to increase significantly Ay (lower 
panel of Fig. |^) for both potentials without affecting appreciably the differential cross-section 
(upper panel of Fig. |^). As the energy increases, it becomes evident that the introduction of 
these 3iVF-like effects provide a good description of Ay in both hemispheres. To obtain this, 
the corrections have to operate differently in the two directions, with a substantial increase 
of the analyzing power in the backward hemisphere, where the peak evolves, and at the same 
time with a slight suppression of Ay in the forward hemisphere, especially for the Bonn B 
potential. Clearly, these irreducible pionic effects have the ability to achieve both goals. In 
Figs. ^ and |^ we compare pd data measured at the proton-equivalent energy of 12 MeV with 
nd calculations at 11.3 MeV (for consistency with the Coulomb slow-down assumption). For 
both potentials, the 3iVF-like effects increase iTu significantly in the backward hemishere 
and suppress slightly the observable at forward angles. Comparison with experimental data 
suggest that these modifications are correct in both directions, although great caution has 
to be exercised when comparing nd calculation with pd data. Around the peak at backward 
angles we observe also a quite large increase for the case of the Bonn B potential, while the 
effect is smaller for the CD-Bonn case. This last feature is similar to what has been observed 
at lower energy. 

The SA^^F-like effects have a very minor impact on T22, also at 12 MeV (upper panel of 
Fig. |]), while it is evident that this observable is more sensitive to the choice of the underlying 
2N potential, and in absence of other additional contributions, we might conclude that the 
data seem to favor more the BBEST calculation, with respect to the CD-BEST results. 

In Fig. P we consider the remaining two tensor analyzing powers, T20 (upper panel) and 
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T21 (lower panel). T20 exhibits an interesting evolution since the results for the BBEST 
+ 0PE-3NF case are significantly different around 110° than the other cases, and they 
are remarkably close to the experimental data. An interesting situation occurs also for 
T21 around 100° where the BEST + 0PE-3NF calculation are appreciably more negative 
than the other calculations considered in the figure. We calculated that this same situation 
evolves also at higher energies (e.g., 18-20 MeV). Unfortunately, no experimental data have 
been found for T21 in this energy range, and new measurements for T21 would be very useful 
in the range of 10-20 MeV. 

As has been explained in Refs. IH], one basic aspect of the irreducible pionic effects 



which generate the 3NF diagram included in our calculation develops as a consequence of 
an "imperfect" cancellation with respect to the mesonic retardation contributions. Since in 
our study we employ realistic nucleon-nucleon potentials, we observe that they are heavily 
based on fitting procedures of experimental nucleon-nucleon data and are therefore a sort 
of "black boxes" with respect to the variety and structure of meson-exchange diagrams 
included, except probably the OPE term, since only the longest range of the nuclear force is 
well established. It should therefore not be a surprise that for each 2A^ potential we included 
in the corresponding 3NF diagram a phenomenological parameter that governs the level of 
cancellation against meson retardation effects. In the approach we developed in Ref. |T3 



we used the experimental value of Ay at the peak to actually fix this parameter. Then, it 
is obviously of interest to study how this parameter evolves with energy. For this reason we 
considered the wealth of pd experimental data for Ay measured at 12, 14, 16, and 18 MeV 
(Ref. |j53|) and compared these data with nd calculations at 11.3, 13.3, 15.3, and 17.3 MeV, 
respectively, finding that we could reproduce how the peak evolves with a perfectly linear 
dependence of the parameter governing the cancellation of the 3NF. 

The situation is shown in Figs. [TO and O for the Bonn B and CD-Bonn potentials. In 
the two figures, the upper panel compares data with results taken without 3NF, while the 
results in the lower panel include 3NF effect. Inclusion of the 0PE-3NF effect provides 
very good results for the Bonn B potential, not only in the region of the peak (around 130° 
MeV), but also in the region around 100° where the dip evolves. For the CD-Bonn potential 
the situation is similarly satisfactorily for the evolution of the Ay peak while the evolution 
of the dip at lower angles is fair but not optimal. For both potentials, we report in Tab. H 
the value of the parameter and the corresponding energy we have employed in calculating 
this 3NF effect. 

In this same energy range, we focus attention on the nucleon-deuteron polarization trans- 
fer coefficients, , , ^ ~yy ^ and ^ . It has been suggested in Ref. that these 
observables are sensitive to the tensor part of the nuclear forces and therefore they could 
in principle represent a good testing ground for the 3A^F-operator we are studying in the 
present paper, since this has a pronounced tensor structure. In addition, recent experi- 
mental data are now available, for both nd and pd systems, and a number of theoretical 
studies about these coefficients have been performed already fS^^Sj with a variety of dif- 
ferent 2A^ potentials, with the addition of 3A^F's, and more recently also with the inclusion 
of the modifications introduced by Coulomb effects. From these studies it emerged that the 
nucleon-to-nucleon transfer coefficients and exhibit a scaling behavior with respect 
to the triton binding energy, while the nucleon-to-deuteron coefficients Ky^''^'^' , and Ky'^' 
do not scale. Moreover, K^' and K^' exhibit sizable Coulomb effects, while for the other 



10 



two coefficients the effects are much less appreciable. For the case of K"^ ■, where it was 
possible to compare directly theory with nd data, it was found that the theoretical calcula- 
tions underpredict the minimum at 110°, once the scaling effect with binding, originated by 
the 27r-exchange 3NF, was properly taken into account. Moreover, discrepancies between 
theoretical calculations and experimental data concerning the Coulomb effects for at 
19 MeV concur to conclude that the situation is not fully understood. For the nucleon-to- 
deuteron tensor transfer coefficients the situation is also unclear, since Coulomb effects and 
traditional 27r-exchange 3NF provide too small modifications for a correct reproduction of 
data (for K^^ the peak at 135° is underestimated). This state of affairs demands for 
more theoretical investigations; at the same time more extensive experimental studies for 
these observables in this energy range could be extremely useful. 

In Figs. |l^ and |13] we show the results obtained with the Bonn B and CD-Bonn po- 
tentials, respectively, and compare these with experimental nd data taken at 15, 17, and 
19 MeV, Ref. |5^. Our calculations do not include the effects of the 27r-exchange 3NF, 



and therefore one should take into account in the discussion the rescaling effect that tend 
to push the lines downward (this tendency is however reduced with the Bonn-types interac- 
tions, which provide smaller 3A^ underbinding with respect to other 2N interactions). The 
results for the Bonn B case suggest that the 3A^F effect we have calculated is able to correct 
the underprediction of the minimum of Kjj at 110°, however the results with the CD-Bonn 
potential do not conffim this indication, since the modifications are smaller here and have 
the tendency to go in the opposite direction. In Figs. |14| and |15| we similarly compare the pd 
data at 19 MeV |5^ with nd calculations at 18.3 MeV. This comparison has to be made with 
great prudence, since Coulomb effects and the 27r-3NF provide appreciable effects, but the 
two modifications tend somewhat to cancel out Nevertheless, it is interesting to observe 
that here the minimum of is appreciably overpredicted in the Bonn B case (lower panel 
of Fig. |T^) while in the case of the CD-Bonn potential the irreducible pionic effects do not 
affect the results appreciably and therefore they maintain the quality of the fit (lower panel 
of Fig. ^). Thus, the experimental "mismatch" between the nd and pd data for at 19 
MeV acquires an interesting twist: the nd data seem to support the Bonn B + 0PE-3NF 
calculations, the pd data seem to support more the CD-Bonn -|- 0PE-3NF results. 



Finally, in Figs. |T6| and O the tensor transfer coefficients are compared with experi- 



ments at 19 MeV. Again, we observe appreciable effects introduced by the 0PE-3NF 
diagram in the case of the Bonn B interaction, while for the CD-Bonn case the effects appear 
to be smaller, although they go in the same direction. We observe that the calculation with 
Bonn B + 0PE-3NF (Fig. 0) is able to solve the underprediction problem in ^ ~yy ^ but 
at the same time it increases the discrepancies in Ky^\ while in the case of the CD-Bonn 
+ 0PE-3NF the situation remains essentially unchanged (Fig. p!7|). 



IV. CONCLUSIONS 

We have studied the low-energy effects in nucleon-deuteron scattering due to an irre- 
ducible pionic effect leading to a 3NF diagram of pronounced tensor structure. Overall we 
find that the effects of this diagram on the differential cross section are negligible in the 
considered energy range, indicating that the average impact of this 3NF diagram on 3A^ 
dynamics is small. On the other hand, for the case of the spin observables, the higher sensi- 
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tivity to the smaller components of the wave function is able to detect the presence of this 
contribution. This is so especially for the case of the vector analyzing powers Ay and iTn, 
which are considered to be a magnifying glass for the triplet P waves of the 2N subsystem. 
Independently of the 2N interactions used as input, we found that this 3NF contribution 
has the potential to significantly increase (about a 30% effect) the magnitude of these two 
observables, and to solve a discrepancy observed long time ago. This is a consequence of 
the specific spin-isospin structure of such 3NF diagram which affects in a privileged manner 
the triplet odd states 

For the remaining spin observables (tensor analyzing powers and polarization transfer 
coefficients) the changes due to this pionic contribution are small, however we found sit- 
uations where the diagram produces appreciable effects. In particular with the use of the 
Bonn B potential there are indications that the corrections produced move towards the right 
direction, but the results with the newer CD-Bonn potential reduce the size of these changes 
considerably. We found few cases where the calculations are slightly - but appreciably - 
different when the 3NF effect is calculated with Bonn B or CD-Bonn potentials. It hap- 
pens with T21 in the energy range 10-20 MeV, with in the energy range 15-19 MeV, and 
with K^^-yy and ' around 19 MeV. There are however still too many uncertainties for 
arriving at a definite conclusion about these slight differences. Comparison of the wealth of 
experimental data for charged particles has been done using only the Coulomb slow-down 
hypothesis, while a consistent inclusion of Coulomb effects is in principle required. Finally, 
a more complete study requires also the inclusion of the remaining 3NF diagrams of dif- 
ferent topology. In particular, it is known that the 27r3NF diagram produces additional 
corrections which are needed for removing the underbinding of the 3N bound state. The 
influence on the vector analyzing powers in the considered energy range, however, appears 
to be marginal. 
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APPENDIX A: BINDING ENERGY 

With the same model interactions used in the main text, we have calculated also the 



triton binding energy. The results are given in the Tab. |T|. For each one of the three 
potentials, the flrst value in the table reports the 3A^ binding energy calculated by a direct 
solution of the homogeneous Faddeev equation in two momentum variables (2D), using the 
original 2N potential as input, and without resorting to the separable expansion method. 
Details on the computational method are explained in Ref. |^ and references therein. The 



second value has been calculated also via a direct solution of the 2D Faddeev equation, 
but this time using as input the separable expansion of the 2N potential, with the same 
ranks as given by Tab. |I[ The third value represents the binding energy obtained with 
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the same separable representation of the 2N potential as in the previous line, but using the 
one-dimensional algorithm corresponding to the Lovelace- Alt-Grassberger-Sandhas homoge- 
neous equation [^. Finally, the last line includes the irreducible pionic effects, as discussed 
in this work, in the ID calculation for the binding energy. 
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TABLES 



NN states 


PEST 


BBEST 


CDBEST 


1 O 

Do 


5 


5 


5(+5) 




6 


6 


6 




5 


4 


4 


'Pi 


5 


4 


4 


'Pi 


5 


4 


4 




5 


4 


4 




5 


4 


4 


'P2 -' F2 


5 


5 


5 



TABLE L Ranks used in the separable expansion for each 2N state and for the various poten- 
tials. For the CDBEST potential in the ^Sq channel, the expansion has been done separately for 
nn and pn cases. 



E (MeV) 


BBEST 


CDBEST 


3.0 


0.730 


0.58 


8.5 


0.733 


0.60 


1L3 


0.743 


0.63 


13.3 


0.753 


0.65 


15.3 


0.763 


0.67 


17.3 


0.773 


0.69 


18.3 


0.778 


0.70 



TABLE IL Energy dependence of the effective parameter used to govern the cancellation in 
Eq. (I). 





Paris 


Bonn B 


CD Bonn 


2D - orig 


-7.385 


-8.101 


-7.958 


2D - EST 


-7.376 


-8.088 


-7.947 


ID - EST 


-7.376 


-8.088 


-7.947 


ID + 3NF 


-7.663 


-7.943 


-8.077 



TABLE in. Results obtained for the triton binding energy (MeV). The first three lines cor- 
respond to different calculational methods without the inclusion of the 2>NF effect. The last line 
includes the effects of the 0PE-3NF diagram. The 2>NF effect has been calculated with the effective 
parameter c determined at 3 MeV. 
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FIGURES 
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45 90 135 



FIG. 1. Differential cross-section and analyzing power for nd scattering at 3 MeV (Lab). Cal- 
culations with the PEST potential (dotted hne), BBEST (dot-dashed) and CD-BEST (dashed). 
For each 2N potential, the corresponding ranks are given in Tab.|. The thick solid line contains 
the resulting modifications introduced by the 0PE-3NF effect, for all three potentials. Data (grey 
squares) from Refs. and Q. 
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FIG. 2. Deuteron analyzing powers T22 (upper panel) and iTn (lower panel) at 3.3 MeV. The 
three thin lines are calculations with 2N potentials only. The corresponding thick lines contain also 
the 0PE-3NF effects. As in the previous figure, dotted, dot-dashed, and dashed lines represent 
respectively the results with PEST, BBEST, and CD-BEST potentials. Data are from Ref. |44|. 
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FIG. 3. Same as in Fig. |2|, but for the T20 and T21 analyzing powers. 
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FIG. 4. Same as in Fig. |^ but for nd scattering at 8.5 MeV. Data are from Refs. |47|, and 
|50|. Calculations are for CD-BEST (dashed line) and BBEST (dash-dotted line) potentials. The 
thick solid line includes the irreducible 3NF-like effects for both potentials. 
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FIG. 5. Same as in Fig. ^, but for nd scattering at 8.3 MeV. Data are for pd scattering at 9 
MeV, from Ref. |^]. Thick (thin) Unes are calculations with (without) inclusion of the irreducible 
pionic effects, as discussed in this work. The calculations are for the CD-BEST (dashed lines) and 
BBEST (dot-dashed lines). 
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FIG. 6. Same as in Fig. |5|, but for the T20 and T21 analyzing powers. 
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FIG. 7. Same as in Fig. ^ but for nucleon scattering at 12.0 MeV. Data for the differential 
cross-section are from Refs. |Q (triangles, nd), (circles, pd), and ||5^ (squares, pd). Data 
for the neutron analyzing power (^4^^) are from Refs. (black circles) and |49| (gray squares). 
Calculations are for CD-BEST (dashed line) and BBEST (dash-dotted line) potentials. The thick 
solid line includes the irreducible 3A^i^-like effects for both potentials. 
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FIG. 8. Same as in Fig. ^, but for nd scattering at 11.3 MeV. Data are for pd scattering at 12 
MeV, from Ref. [^]. Thick (thin) hnes are calculations with (without) inclusion of the irreducible 
pionic effects, as discussed in this work. The calculations are for the CD-BEST (dashed lines) and 
BBEST (dot-dashed lines). 
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FIG. 10. Evolution of the analyzing power Ay in the range 12-18 MeV. Squares, circles, dia- 
monds and triangles represents pd data at 18, 16, 14 and 12 MeV, respectively, taken from Ref. |53|. 
The lines in the upper panel refer to corresponding calculations with the BBEST potential, while 
in the lower panel the 0PE-3NF effects are also included. 
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FIG. 12. Spin-transfer coefficient -fC|' at 15 MeV (upper panel), 17 MeV (middle panel), and 
19 MeV (lower panel). Comparison between nd data from Ref. and BBEST calculations. Solid 
(dashed) lines include (exclude) the irreducible 3A^F-like pionic effects. 
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FIG. 13. Same as in Fig. 12 but for the CD-BEST potential. 
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FIG. 14. Calculations for the nd spin-transfer-coefficients (upper panel) and K'jj' (lower 
panel), at 18.3 MeV, for the BBEST potential. Data are for the equivalent pd observables at 19.0 
MeV, from Ref. |54]. Solid (dashed) lines include (exclude) the irreducible 3A^F-like pionic effects. 
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FIG. 16. 



Polarization transfer coefficients Ky ^ 



-y y 



(upper panel) and Ky ^ (lower panel) 
calculated at 18.3 MeV for the BBEST potential. Data are for pd scattering at 19.0 MeV, from 
Ref. |58|. 
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